Far infrared giant dipole resonances in neutral quantum dots 
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Abstract 



A resonance behaviour of the far infrared absorption probability at a fre- 
quency ~ iV 1 / 4 is predicted for clusters of N electron- hole pairs (2 < N < 110) 
confined in disk-shaped quantum dots. For radially symmetric dots, the ab- 
sorption is dominated by a Giant Dipole Resonance, which accounts for more 
than 98 % of the energy-weighted photoabsorption sum rule. 
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I. INTRODUCTION 



A well-established result of Nuclear Physics is the presence of giant resonance modes 
related to collective excitations in nuclei. Particularly, collective states corresponding to 
giant electric dipole oscillations are known in nuclei long agc0i. They show themselves as 
very high peaks in photoabsorption at energies ~ SOA" 1 / 3 MeV, where A is the mass number. 
An understanding of this phenomenon can be achieved by means of a simple hydrodynamic 
model in which one considers small density perturbations against the background nuclear 
density and computes the small-oscillation frequencies associated to dipole excitations!. A 
second, microscopic, approach makes use of the nuclear shell model@. It is shown that the 
giant dipole resonance (GDR) is a state coming from the splitting of a degenerate set of 1~ 
states, which takes almost the whole strength of the 1~ transition. 

Collective states very similar to the GDR have been extensively studied, both exper- 
imentally and theoretically, in Atomic Physics. They are at the base of giant and broad 
atomic resonances in the photoionisation continuum!. The photoabsorption cross section of 
metallic clusterslH exhibits also the dominance of a collective oscillation mode at 2 - 3 eVi. 
Theoretical calculations based on the Random Phase Approximation (RPA)!@ are shown to 
reproduce well the resonance position and dipole strength distribution. 

In the present paper we show that, in electron- hole systems confined in a quantum dot, 
the energy-weighted photoabsorption sum rule is almost saturated by GDR-like states, in 
which the electronic cloud oscillates relatively to the hole cloud. This result is mainly related 
to the existence of positively and negatively charged particles, and not to a particular dot 
geometry. Thus, we will use a very simplified model of disk-shaped dot, in which the lateral 
confinement is ideally parabolic!. 

N electron- hole pairs are supposed to be created in the dot by means of, e. g., laser 
pumping with a typical ~ 2 eV energy, as corresponding to a semiconductor band gap. The 
created particles live hundreds of picoseconds and even more0. We will study the linear 
absorption of a second far infrared (FIR), ~ 5 meV, light wave by the A^-pair cluster. The 
period of this signal is ~ lCT 13 s, thus the cluster may be considered as stable when studying 
absorption. We will use the dipole approximation for the interaction hamiltonian. 

Experiments on FIR study of excitonic systems have been conducted for a long time0. 
More recently, the modifications of the interband photoluminescence (PL) of quantum well 
excitons caused by the absorption of FIR radiation have been studiedtj by means of a 
technique baptised as ODR (optical detection of FIR resonances). These experiments are 
conducted at relatively small exciton densities, where many-body effects are not expected 
to play a significant role. To study the collective FIR absorption by a multi-exciton system, 
intense laser beams for the interband processes shall be used. Experimentally, densities as 
high as 10 12 electron-hole (e-h) pairs/cm 2 in GaAs wells have been achieved by using pulsed 
(5 ps long) dye lasers^!. The effective absorbed power is greater than 2 MW/cm 2 . 

To our knowledge, however, the effect of FIR radiation on dense excitonic systems have 
not been experimentally studied so far, although intense enough laser beams and the ODR 
technique are available. For simplicity of the calculations, we confine the system in a quan- 
tum dot (qdot). This confinement makes the detection a little more cumbersome, but still 
experimentally achievable. Single qdot multi-excitonic PL has been isolated from the aver- 
aged PL by means of microscopy techniques (see, for example, Ref. [|l4jl). Thus, in principle, 
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the ODR study of single qdots is also possible. 

The organisation of the paper is as follows. In Section II, we present the model to be 
used in the paper. Next, the basic RPA equations are explicitly written in Section III. In 
the later two sections, the main results and a few concluding remarks are given. 

II. THE MODEL 

We will study a quasi-planar quantum dot, i. e. very narrow in the growth direction (the 
z axis), in such a way that the first quantum well sub-band approximation works, and no 
significant modifications of Coulomb interactions due to z-averaging arise. A direct band- 
gap material and high barriers for both electrons and holes in the z direction are assumed 
to exist. Confinement makes also the hole sub-bands to split. Thus, only one e and one h 
bands will be considered, e-h exchange will be neglected. As mentioned above, the lateral 
confinement will be modelled by an ideal parabolic potential. Finally, for simplicity, we take 
m e = rrih. This is a non realistic assumption, justified only by the absence of experimental 
results. In a GaAs dot, for example, 8 nm wide in the z direction, the ratio of the in plane 
heavy hole to electron masses is around 2. The equal- mass simplification makes calculations 
more easy, giving nevertheless a qualitatively correct picture of what is happening. 

The second-quantised Hamiltonian for the two-dimensional motion of electrons and holes 
is written in oscillator units as 

H = ^ej(ejei + h]hi) + 7^ E (hj\yr\k,l)4e}eie k 

+f E Ml/r|M)^]U fc -/3 E (i,j|l/r|M)^Afe, (1) 

i,j,k,l 

where ei(h%) are electron (hole) operators, (3 = J , k is the dielectric constant, 
and u is the dot frequency. The Coulomb matrix elements are defined as 

— ; ^0K^l)0K^2)0fc(ri)<M^2), (2) 

\T\ — r 2 \ 

and the explicit form of the harmonic-oscillator orbitals is the following 

h = C nj!mj r^L^(r 2 )e- r2 / 2 e im ^ e , (3) 

where rij and rrij are oscillator quantum numbers, and C n> | m | = Jn\/[ir (n + |m|)!]. The 
one-particle energies in the quadratic potential are 6j — 1 + 2rij + \rrij\. 

We consider the dot in the strong-confinement regime, where the one-particle energy 
dominates the total energy. It means, (3 < /3 C , where (3 C may be estimated^ from the 
breakdown of perturbation theory or from the onset of BCS pairing to be (3 C w 0.58 N 1 / 4 . 
The associated value of density is p c ~ 4/ (ira 2 B ). For GaAs, the Bohr radius is as ~ 14 nm, 
thus p c ~ 5 • 10 11 pairs/cm 2 . 

For simplicity, we study closed-shell quantum dots. The number of electrons is equal 
to the number of holes, and is given by the expression iV = N S (N S + 1), where N s is the 
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last harmonic-oscillator filled shell. The zeroth-order ground state-function is a product of 
Slater determinants for electrons and holes. The total angular momentum, total electron 
spin and total hole spin are all equal to zero for this state. 

The lowest excitations against the ground state are one-particle excitations from the last 
filled shell to the next empty shell. The states with total angular momentum, M — ±1, and 
electron and hole total spins, S e = Sh = 0, are the candidates among which we shall look 
for the GDR. To be definite, we will consider M — 1. Notice that there are N s orbitals in 
the last shell, that is 2N S electrons and, thus, 2N S wave functions with M = 1. There are 
also 2N S functions corresponding to hole excitations. It makes a total of 4N S many-particle 
states. 

The Coulomb interaction breaks the oscillator degeneracy. The collective nature of the 
GDR is manifested in two facts: (i) it acquires the highest excitation energy in the multiplet, 
and (ii) the electric dipole transition probability from the ground state is strongly enhanced!. 
We will present RPA calculations in support of these statements. 



III. THE RANDOM PHASE APPROXIMATION 

In the RPA0, we allow a general correlated ground state, \RPA), and the excited states 
are looked for in the form 

V = Q^\RPA), (4) 
where the operator is given by the expression 

Q ] = E(^a4ca + X} x hlh x - Y* x e{e a - Y*AK\ (5) 

(7, A 

The index A runs over the last Hartree-Fock (HF) filled shell, and a runs over the first 
unfilled shell. The X coefficients are nonzero only for "up" transitions, i.e. for which the 
angular momentum values satisfy m a — nix — 1, and the Y coefficients are nonzero only for 
"down" transitions, i.e. A and a satisfying ni\ — m a = 1. They are obtained from the RPA 
equations 

A B \ fX\ A „ (X 
-B l 

in which the A matrix is given by 



B < -A')\yI- AE (yI- < 6 > 



/ A ee A eh 
\ A^ e 



(7) 



and its matrix elements are the following 



K%r, = - ff)t^ +P((cr, l*\l/r\\, t) - (a, f,\l/r\T, A)) , (8) 
Af x ^ = -(5{a^\l/r\\r). (9) 

A = A ee , A he = A eh because of the electron- hole symmetry of our hamiltonian. r and 
H are indexes similar to a and A, respectively. The B matrix has a similar structure. Its 
matrix elements, calculated between up and down transitions, are given by 
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B?x, T , = /3((a,r|l/r|A,/i) - (cr,r|l/r|//, A)), (10) 



Sg )T ^ = -/3((T J T|l/r|A l /i>. (11) 

^h/i _ gee an( j ^fi,e _ ^e/i j g transpose of the numerical matrix B, and the A d 
matrix is calculated from (|S],|9]) between down transitions. Usually, positive (physical) and 
negative (unphysical) excitation energies of equal magnitude come from (||). The physical 
solutions annihilate the RPA ground state 

Q\RPA) = 0, (12) 

and satisfy the normalisation condition 

1=E{(^a) 2 -(F ctA ) 2 }. (13) 

<7.A 

The FIR radiation to be absorbed is assumed to be linearly polarised, with the electric 
field vibrating along the y direction. The dipole matrix elements are computed from 

(*\D\RPA) = £{PCa - X* x )(a\y\\) + (Y a \ - Y a \)(X\y\a)}. (14) 

(T.A 



The calculated dipole elements fulfil the energy-weighted sum ru 

A£ mD\RPA)\ 2 = £ AEoK^ol^lO)! 2 = N/2, (15) 

where indexes refer to the noninteracting f3 = case. 

The HF energies and one-particle wave functions were obtained selfconsistently by means 
of an iterative scheme. We write 

\a) = J £C ai \i), (16) 

where latin indexes refer to harmonic oscillator orbitals. 272 orbitals, i. e. 16 harmonic 
oscillator shells where included in the calculations. The HF equations read 

J2i e i 5 v J2 C ^ C Ld(i,k\±/r\l,j)}C aj = e a C ai , (17) 

j n<F k,l 

and the starting approximation is C a i = 5 a i- F denotes the Fermi level. A few iterations are 
enough to reach convergence at very low (3 values, while at f3 ~ 1, 15 iterations are needed. 



IV. RESULTS 

Excitation energies and dipole strength distributions were computed for N ranging be- 
tween 2 and 110, and (3 in the interval < (3 < 0.6 iV 1 / 4 . For larger values of /?, pairing 
becomes very importantll!, and should be taken into account more properly. Let us stress, 
however, that pairing should reinforce the picture of GDR dominance of photoabsorption, 
leading to an increase of the excitation energy and, consequently, to a decrease in the same 
proportion of the dipole element squared. 
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The results for N = 110 pairs at (3 = 1 are shown in Fig. [I], where dipole matrix elements 
squared vs excitation energies are shown. The fraction of the sum rule flTBD accounted for 
the GDR state is close to 100 %. Notice that only collective states, for which X e = —X h , 
are represented in the figure. The rest of the states show a threefold degeneracy and give 
practically zero dipole matrix elements. Let us note that the energy gap is increased (AE > 
1) after the Coulomb interaction is switched on. This is a signal that the overall interaction 
among particles is attractive. 

We draw in Fig. Qa the GDR excitation energies at (3 = 1 as a function of iV 1//4 . Notice 
the approximate dependence AE gdr — 1 « 0.7 + 0.4 iV 1//4 . At f3 ^ 1, the dependence is 
AEgfo — 1 ~ (0.7 + 0.4 N l / 4 )(3. Turning back to ordinary units, we get an excitation energy 



AE gdr ~ (1 + 0.7 p)hu + 0.4 J2tt/3 p 1/2 , (18) 

mas 

where p ~ ?>N l l 2 moj / '(2irh) is the density of pairs. For large enough N, and u — » in 
order to preserve a constant density, AE gdr m 0.6 h 2 p 1 / 2 / '(mag), leading for example to 
AE g( ir ~ 4.5 meV for typical values of the constants in GaAs and p ~ 5 ■ 10 11 pairs/cm 2 . 

The GDR dipole elements squared at (3 = 1 as a function of N are shown in Fig. gb. 
The points fit a linear dependence ~ 0.17 N. For (3^1, the saturation of the sum rule 
implies that D 2 « N/{2AE gdr ) « iV/(2 + (1.4 + O.SiV 1 / 4 )/?). 

Finally, we draw in Fig. |3] the fraction of the energy- weighted sum rule flTSD, the 51 sum 
rule, corresponding to the transition to the GDR state. This fraction is higher than 98% 
for any of the systems shown in Figs. |3|a and ^|b. (We got a slight violation, fraction=1.01, 
for iV = 110 at (3 — 1, but this is an artifact related to the finiteness of the basis used to 
construct the HF states). 



V. CONCLUSIONS 

We have shown that the far infrared absorption spectrum of neutral confined systems of 
electrons and holes is dominated by a giant dipole resonance which accounts for more than 
98% of the energy-weighted photoabsorption sum rule. Calculations were done in a very 
simplified two-band model with m e = rrih and a disk-shaped parabolic dot. The qualitative 
conclusions are, however, expected to be valid for realistic systems since they are mainly 
related to the existence of positively and negatively charged particles in the system, which 
causes the enhancement of dipole oscillations. 

For the resonance position at intermediate N, we obtained the fit, AE g d r = {1 + 
(0.7 + 0.4 N l / 4 )f3}Tiuj. In terms of the density, we obtain AE gdr w (1 + 0.7(3)huj + 

0.4^271-/3 Ti 2 p l / 2 /(m os). The dipole strength scales with N: D 2 ~ N, at small N val- 
ues, and D 2 ~ iV 3//4 at larger values of N. 

The GDR state we found is the analogue in a finite two-dimensional system of the out- 
of-phase plasma frequency in the two-component plasma0. The term proportional to htu in 
the excitation energy is a shell effect of the external harmonic-oscillator potential. 

We worked in the strong-confinement regime, (3 < 0.6 iV 1/4 . For still higher values of j3, 
BCS pairing should be explicitly included in the calculations. Recently, we have computed 
the ground-state energy of the confined iV-pair system by using a BCS trial functionEj. The 
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effect of pairing in AE g d r could be taken into account by means of a quasiparticle RPA 
schemeB, which may be taken as the analogue for a finite system of the Gor'kov's pairing 
approximation in the excitonic insulator phaseliHil. At relatively high densities we expect, 
qualitatively, a shift of the resonance position in accordance to the opening of the gap at 
the Fermi level: 



AE gdr = \ {2A F ) 2 + ( (1 + 0.7P)tiu + 0.4^2^/3 Ti 2 p l l 2 /(m a B )) , (19) 



where is the gap function at the Fermi energy The dipole strength shall decrease in the 
same proportion, as dictated by (|T5|). 
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FIGURES 

FIG. 1. Dipole matrix elements squared for 110 pairs at (3 = 1. 

FIG. 2. a) GDR excitation energy as a function of the number of pairs in the dot at (3 = 1. 
(b) GDR dipole matrix elements squared as a function of N at (3 = 1. 

FIG. 3. Fraction of the energy- weighted sum rule accounted for the GDR: a) As a function of 
(3 for iV = 110. b) As a function of the number of pairs at (3 = 1. 
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